1 



ORNL-CTP-9703 and lhep-ph/9712320 



Heavy Quarkonium Production and Propagation in Nuclei 
Cheuk-Yin Wong^* 

^ ■ ^Physics Division, Oak Ridge National Laboratory, Oak Ridge, TN 37831 
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^ ■ We describe a precursor in heavy quarkonium production in terms of a coherent ad- 

D ■ mixture of states of different color, spin, and angular momentum quantum numbers, and 

Q . obtain the production amplitudes for different quarkonium bound states by projecting 

Q^ ! out this precursor state onto these bound states. The precursor is absorbed in its passage 

^ I through a nucleus in a pA reaction, and the total cross section between this precursor with 

^ ■ a nucleon can be calculated with the two-gluon model of the Pomeron. Such a description 

O . of coherent precursors and their subsequent interactions with nucleons can explain many 

^ . salient features of J/ ip and ip' production in pA collisions. 

(N : 

^ , 

■ 1. Introduction 

as ; 

^ ■ It is a great pleasure to write this article in honor of Prof. Tai-You Wu's 90th Birth- 

Ph, day. Prof. Wu's grandfather and father were respectively "Jin Shi" and "Ju Ren" , which 

Oh! were among the highest honors a Chinese scholar could achieve in the last century. 

^ I Prof. Wu himself was instrumental in bringing modern physics to China, and has trained 

L| ■ a large number of physicists, including C. N. Yang and T. D. Lee. Through the lineage 



X 



of Tai-You Wu, Chinese physics scholars of today are a part of the continuation of the 
Chinese scholarship of past few thousand years, and the solid record of accomplishments 
5^ I of past Chinese scholars will inspire Chinese physicists to continue to make contributions 

to the world of science in the future. 

In the search of for the quark-gluon plasma, it has been suggested that the production 
of charmonium will be suppressed in a quark-gluon plasma because of the screening of 
the interaction between c and c [0] . To extract information on the suppression due to the 
quark-gluon plasma, it is necessary to study the suppression of J/ ip production by sources 
different from the quark-gluon plasma. It is therefore useful to examine the mechanism 
of heavy quarkonium production and its propagation in nuclei [pHTl]] . 

To follow the processes of production and propagation of a charmonium or a heavy 
quarkonium, one needs to understand the nature of the object that is formed by the 
collision of partons in a hard-scattering process. This leads to the picture of a precursor 
formed by parton collisions that is a coherent mixture of different color, spin, and angular 



*This research was supported by the Division of Nuclear Physics, U.S. D.O.E. under Contract DE-AC05- 
960R22464 managed by Lockheed Martin Energy Research Corp. 



2 



momentum states. The production and the subsequent interaction of this object with 
nucleons along its passage through a nucleus is the subject of the present study. 

The present work is the result of a continuing evolution and synthesis of many previous 
attempts to understand production and the absorption of heavy quarkonium in a nucleus. 
Previous studies on heavy quarkonium production include the color evaporation model 
[p!^ , the color-singlet model |jl3l, and the color-octet model fl^. The absorption of 
J/ip has been previously described as the propagation of a hybrid object [(00)8(7] 0) the 
additive quark model , a color-singlet coherent state due to the multiple collisions with 
nucleons ||,^ , and an incoherent admixture of color singlet and color-octet states [p!0| , p^ . 
The model we shall present in terms of a coherent color admixture of the precursor provides 
a description of the important features of the production and absorption process. There 
are also specific predictions of this model which can be further tested by experiment. 



2. Production Process 



We can study the collision of the parton 6 of a beam nucleon with the parton a of the 
target nucleon leading to the production of a Q-Q pair which later materializes directly 
or indirectly into a bound heavy quarkonium. The initial state of the Q-Q pair from 
the collision at tj = is represented by the state vector 

\^ab{ti) >= M{ab ^ Q{P/2 + q)Q{P/2 - q))\Q{P/2 + q)Q{P/2 - q) > (1) 

where J^{ab Q{P/2 + q)Q{P/2 — q)) is the Feynman amplitude for the a + b ^ Q + Q 
process, P is the center-of-mass momentum, q is the relative momentum of the Q-Q pair. 
For simplicity of notation, the color and azimuthal spin components are understood. 
One can perform a decomposition in terms of color and angular momentum states as 

l^abiU) >= E ^SLsiQ)\QQfLS]iP) > (2) 

CJLS 

where 

0JL5(?) =< QQl'LSmiMiab ^ QiP/2 + q)Q{P/2 - q))\Q{P/2 + q)Q{P/2 - q) > (3) 

with \QQ[^ L'j]{P) > describing the center-of-mass motion of the QQ pair in color state 
C, and angular momentum quantum numbers JLS. 

The state will evolve according to perturbative QCD, 

I'^abit) > = U{t,U)\<l>ab{ti) > = T exp{-t f Hidt} \^ab{ti) > (4) 

Jti 

where the evolution operator U(t,ti) can be expanded out in a time-ordered perturbation 
series in terms of the interaction Hj . The state $a6 can be called the heavy quarkonium 
precursor state formed by partons a and b. Because the partons are propagators whose 
legs joined onto their parent nucleons, the partons are off the mass shell, with energies 
that depend on the invariant masses of the other products of the nucleon-nucleon collision 
(see [0 and Eq. (4.14) of ||16|). Thus, the energy of the state ^ab also depends on the 
invariant masses of the other products of the nucleon-nucleon collision. It is reasonable 
to associate the on-shell energies of the state ^ab as an average energy of an admixture 
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of bound states of different energies and obtain the bound state production amplitude by 
projecting the asymptotic precursor state at t = oo onto various bound states within a 
certain energy range. For example, guided by the evaporation model, one can consider 
this range as between 2mc and 2mD for cc quarkonium production. 

We describe a bound QQ state with the quantum numbers JLS and other quantum 
numbers with a center-of-mass momentum P as 



\'^jLs;Pq >-- 



^^'''-RjLs{q)\QQfL^J^m>. (5) 



\ ArnqTriQ' 



The above state may be specified by additional quantum numbers such as the number of 
radial nodes, etc. For simplicity of notation, these additional quantum numbers will not 
be written out explicitly. 

In the lowest-order perturbative QCD, the probability amplitude for the direct produc- 
tion of JLS is obtained from projecting Tahiti) onto ^ jls- The projection is simplest in 
the Q-Q center-of-mass system where P = {Mjls, 0) and q = (0, q) and the probability 
amplitude is [0,0 



< ^JLS; Pql^abit = OO) > lowest order = < JLSl Pq\^abiti) > 

73 



\j imgrnQ J (27r)' 

where Mjls is the mass of the bound state, rriQ the mass of the quark, and we have 
followed the normalization of Ijl^. Because the bound state "^jls is a color-singlet state, 
the above projection will involve only color-singlet components of the admixture in Eq. 

®- 

In the next-order perturbation theory, the color-singlet bound state "^jls accompanied 
by a gluon g can be produced by the color-octet component of $a6 in Eq. (^. The 
probability amplitude for the production of the bound state "^{JLS) accompanied by a 
gluon is 

< [^jLs;Pq]9\^ab{t) >=< [^jLs;Pq]9\U{t,ti)\<!>ab{ti) > . (7) 

The state "^jls can also be produced indirectly through the production of different 
bound states "^j'us' which subsequently decay into "^jls- For example, in J/ip produc- 
tion, a large fraction of the observed J/tp comes from the radiative decay of xi and X2 
[0. The production probability, including direct, indirect, and color octet contributions, 
is then the sum of the absolute squares of various amplitudes. 

Heavy quarkonia can be produced by different parton combinations such as g-g, 
q-q, and g-q collisions, which will lead to different precursor states. The total production 
probability will be the sum from all different precursor states. 

For lack of a better name, we shall call the present model the coherent precursor model. 
Such a model has many features similar to the color evaporation model. In fact, one can 
interpret the phenomenological fractional coefficient in the color evaporation model as the 
absolute square of various projection amplitudes. The approximate success of the color 



evaporation model [E^ also assures that this projection treatment will lead to a good 



description of many pieces of experimental data. 
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3. Interaction of a coherent precursor with a nucleon 

The interaction of the precursor with hadron matter or other medium depends on the 
relative kinetic energy of this precursor with respect to the hadron matter or the medium. 
For processes when the precursor travels at a high energy relative to hadron matter, the 
interaction of the precursor with respect to a hadron in the medium can be described 
in terms of Pomeron and Reggeon exchanges. On the other hand, at energies where the 
precursor is nearly at rest with respect to the medium, the interaction depends on the 
binding energies of various bound states. 

We shall consider here the case of charmonium production in pA collisions at a fixed 
target energy of a few hundred GeV fUJH, and T production at 800 GeV At 



these energies, the production of a heavy quarkonium state with xp > will involve the 
interaction of the precursor with nuclear matter at high energies. The total cross section 
between the precursor and a model meson M{q2,qi) can be obtained by using the two- 
gluon model of the Pomeron p6|-^,p|] . Using the notation of p8|, the elastic scattering 



amplitude for the exchange of two gluons between the precursor ^ab{QiQ2) and the meson 
M{q^q4^) in this model is 

Ais,t) = '-j^l f/be-^Q'^ < ^abiQiQ2)M{q^q,)\V\^abiQiQ2)M{qsq,) > (8) 

where s = (^($^6, initial) +p(M(g3g4), initial))^ t = final) -^($^6, initial))^ = Q^. 

For the exchange of two-gluons, V is p8[ 

V = E[A?A^^i3 + i\^2f{XirV,, - XUKfVu - i\^2fXtVi3?, (9) 

a 

where A" is the generator of SU(3) associated with the ith quark, Vij is the interaction 
between the i and the jth quark or antiquark. Because (^3^4) is coupled to a color-singlet 
state, the matrix element in the above equation between different color states vanishes, 
and we have 

A{s,t) (10) 

where -R(r34) is the 53-^4 relative wave function in the model meson M{qsq4). One 
can use the optical theorem to obtain the total cross section for this coherent pre- 
cursor. One can compare the total cross section '7j£,5(o6) obtained from the forward 
elastic scattering amplitude for the case if $ab is a (properly normalized) pure state of 
^jLsi^i'i)\iQiQ'i[^^j]iR) >• Then Eq. ([10|) gives the total cross section between the 
precursor (formed by the partons a and b) given by 

atotiab) = Y: f?Lsi^b) aSr^siab) (11) 

CJLS 

where = / dri2\^SLsi^i2)\V < ^ab\^ab > • (12) 

The above describes the cross section for a coherent precursor colliding with a meson. 
The result of Eq. (O) can also be generalized to that for the collision between the 



precursor and a nucleon as is done in [28 
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From the work of Dolejsi and Hiifner pS), and Wong [^j, we know that the (QQ)^^^- 
nucleon total cross section for a pure color-singlet {QQY^^ state is approximately pro- 
portional to the root-mean-square of the separation between Q and Q. The coherent 
precursor produced after the hard-scattering process has a small spatial extension, so it is 
reasonable to take the color-singlet cross sections a^H for different JLS states to be small 
and nearly equal. The {QQY^^ -imcleon total cross section for a pure color-octet (QQY^^ 
state is insensitive to the Q-Q separation. This approximate independence of the color- 
octet cross section with respect to the Q-Q separation |2^J^ allows us to approximate the 

(8) 

color-octet cross sections for different color-octet JLS states by a single aiJ. Therefore, 
the precursor-nucleon total cross section is 

^totiab) ^ f'\ah) a« + f'\ah) aLI, (13) 



where f {ah) = ^ fjLsi^^b). (14) 

JLS 

In heavy quarkonium production, there can be different combinations of partons leading 
to different precursor states, and each precursor has a total precursor-nucleon cross section 
given by an expression such as Eq. (plsl). Thus, the precursor cross section depends on the 
color-octet fraction of the precursor. 

The most important parton collisions for heavy quarkonium production are the gg and 
qq collisions. In gg collisions, the coherent precursor can be in color-singlet and color-octet 
states, while in qq collisions, the precursor is a color-octet state in the lowest order. For 
the total yield at fixed target energies, the dominant production process comes from the 



collision of gluons |22 



In the kinematic region where the collision of the precursor with nucleons takes place 
at high energies, as in the production of heavy quarkonium in the region of xp > 
in pA collisions, the binding energy of the heavy quarkonium is much smaller than the 
precursor-nucleon collision energy. The probability for the breakup of the Q-Q pair into 
Qq and qQ by combining with light quarks is large in a precursor-nucleon collision. It is 
reasonable to assume that in high-energy precursor-nucleon collisions the absorption (i.e. 
breakup) cross section, aabs{cib), is approximately given by the precursor-nucleon total 
cross section atot{,ob) of Eq. (|T3]). 

4. Propagation of a coherent precursor in nuclear matter 

We consider a pA collision in which a heavy quarkonium precursor is produced in one 
of the nucleon-nucleon collisions and the precursor subsequently travels in the nuclear 
medium of the target nucleus. Because of the coherence of the precursor, it propagates 
through the nuclear medium as a single object with a single absorption cross section 
after its production. When the precursor travels inside the nuclear medium, the time of 
evolution t can be represented equivalently by the corresponding path length L/w, where 
V is the velocity of the precursor in the medium. The state vector after propagating a 
distance L in the nuclear medium is related to the state vector after production by 

|<l>„5(^) >= e-^'^-»("'')^/2|^,,(L = 0) > , (15) 
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where aabsio-b) = crabs{^ab—N) is the precursor absorption cross section for the coUision 
of the precursor with a nucleon, and p is the nuclear matter number density. 

After the passage through a path length of L in the nuclear matter, the amplitude for 
the production of a bound state is obtained by projecting the state vector ^ab{L) 

onto the bound state "^jls- Therefore, 

< ^jLs; Pq\<^ab{L) >= e-P-^^^^-^)^!-^ < s^j^^- Pq\^,,{L = 0) > (16) 
and for the octet production of \^ jls] P<1 > with the emission of a gluon 

< [^jLs;Pq]g\U{t,L/v)\<!>ab{L) >= e-'"^-^^''^)^/^ < j^s; Pq]g\U{t,0)\<^ab{L = 0) > 

When we include precursors from different parton collisions leading to the production of 
the bound state "^jls, there will be different survival factors e~^°'°-'"'^"'^^^^'^ for different 
parton combinations a-b. At fixed-target energies, where the total yield of J/tp in the 
forward direction is dominanted by contributions from gg collisions ||22|, one expects that 



there is essentially only a single survival factor for the total yield in forward directions. 

5. Consequences of the coherent precursor model 

One can outline a few distinct consequences of the model we have presented. In the 
situation when only a single absorption combination of partons dominates the production 
process, the ratio of the production of various bound states in a pA collision will be 
independent of the mass number of the nucleus. This arises because the precursor is 
absorbed in its passage through nuclear matter by a single precursor-nucleon cross section 
o'absiO'b), and the production of all different bound states comes from the projection of the 
precursor state onto the bound states after the absorption. For example, as a consequence 
of Eqs. dll). 



<^'|$33(L) > |2 _ I < ^'|$,,(L = 0) > 



<^|$,3(L)>|2 I < = 0) > 



:i8) 



There is a similar relation for the production amplitude from the color-octet amplitudes 
because of Eq. (p!?!). Therefore, we have 

c^miL) ^ a(V;0(L = O) ajpA ^ jj' X) ^ ajpp ^ j;' X) 

(j{il)){L) a(V^)(L = 0) °^ a{pA^i)X) a{pp i) X)' ^ ' 

which is independent of A. A similar independence with A is predicted for T'/T. Exper- 
imentally, the ratio ip' /{J/il)) is observed to be approximately a constant of the atomic 
numbers ||2^, in agreement with the present picture. The experimental yields of T and 



T' have considerable uncertainties, but within the large experimental uncertainties the 
ratio a{pA — > T' X)/(j{pA T X) is approximately independent of A |^|. The present 



model predicts further that the ratio of the x yield to the J / ijj yield will also be indepen- 
dent of the mass number in pA collisions. It will be of interest to test such a prediction 
in the future. 

A similar model of coherent state production was put forth earlier by Dolejsi and Hiifner 
and Wong [Q, with some major differences. In these models, a color-singlet coherent 
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state of J /if) and ip' is assumed to come from the successive multiple scattering of the 
precursor with a row of nucleons, whereas in the present model, a coherent state of both 
color, angular momentum, and other quantum numbers is formed in the hard-scattering 
process before multiple collisions with nucleons and collisions with nucleons at high energy 
are assumed to lead to absorption, and not to the change of the coherent mixture. 

We note that the total cross section atot for a coherent color admixture is the weighted 
sum of the corresponding cross sections for pure color-singlet and color-octet states. For 
a coherent state that is formed after the hard scattering, the spatial extension of the 
state is small. One expects that o"£] should be nearly zero while 0"^], which is nearly 
independent of the radial size of the state, is very large, of the order of 30 to 60 mb. 
Thus, an absorption cross section of the order of 6 mb extracted from the experimental 
survival factor in pA data 0,^,0 for J/ip and ip' can be understood as arising from a 
color-octet fraction /^^^^ of the order of 10% to 20%, as given by Eq. (|T^. 

Experimentally, the production cross sections for xii X21 and direct J/ip are 131 nb, 
188 nb, and 102 nb, respectively, for vr-N collisions at 300 GeV [^]. Therefore, direct J /il^ 
production is about 25% of the total yield of charmonium. As Xi and X2 are produced 
predominantly via the color-singlet states, while the direct J / if) production comes predom- 
inantly from color-octet states pO| , |2T|] , it is therefore encouraging that the fraction of 25% 
of direct J/ip among the total charmonium yield is of the same order as the color-octet 
fraction /*^^) of about 10%-20% that is consistent with the absorption cross section and 
the theoretical estimates from the two-gluon model of the Pomeron. It will be of interest 
to calculate this color-octet fraction Z*^^-' from Eqs. (^) and (|Tip theoretically. This may 
rectify a previously unsatisfactory outcome when the precursor is described as a mixed 
state with an incoherent color admixture. The incoherent color admixture description of 
the precursor is not satisfactory because the absorption cross section for the color-octet 
state is found, in this picture, to be about 15 to 20 mb in [|l^ and 11 mb in [|ll|], too small 
to be consistent with the theoretical estimates of about 30-60 mb from the two-gluon 
model of the Pomeron p8|,|9|. 



Finally, in different kinematic regions, the importance of the different parton production 
components will change [Q. For example, while gg fusion dominates in the region of small 
xp, which gives the largest yield, the yield at xp » 0.6 is dominated by q-q collisions. 
Unlike the gg fusion which can lead to both color-singlet and color-octet components, qq 
collisions lead to color-octet in the lowest-order approximation. The present model will 
predict an increase in the absorption cross section in the region of large xp due to the large 
color-octet precursor-nucleon cross section. Experimental data at different xp show an 
increase in the effective absorption cross section as xp increases [p!0| , p3| , p^ . However, the 
effect of energy loss needs to be taken into account at these kinematic regions to confirm 
the effect that the color-octet state leads to a very large increase in the absorption of the 
precursor at xp > 0.6. 

In the transition region when two different precursors contribute about equally, as in the 
region around x/ ~ 0.6, the survival probability is determined by two different exponential 
survival factors, with two different absorption cross sections for the two kinds of precursors, 
as a function of the nuclear matter path length. Future experimental investigations on 
the change of the characteristics of the mass dependence from small Xf to large xj will 
be of interest. 
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In conclusion, the coherent precursor model provides a reasonable description of the 
production and propagation of the object form before the formation of the bound states. 
Further tests of the model can be carried out to confront with experiment. 

The author would like to thank Dr. C. W. Wong for comments and communications. 
He also wishes to thank Drs. R. Blankenbecler, H. W. Crater, J. C. Peng, and J. S. Wu 
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